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From the particle physics point of view, the most peculiar property of the dark matter particle is its 
stability on cosmological time scales. We briefly review the possible origins of this characteristic 
feature for candidates whose relic density results from the thermal freeze-out of their annihilation. 
We emphasize that each stabilization mechanism implies an all specific phenomenology. The 
models reviewed include supersymmetric and non-supersymmetric models where the stability is 
a consequence of grand-unification, models where stability is due to an unbroken gauge group 
and models where the DM stability is accidental. The latter possibility includes minimal dark 
matter, hidden vector dark matter and composite DM models. 
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1. Introduction 

The dark matter (DM) particle is fantastically stable. Its decay lifetime has to be obviously 
larger than the age of the Universe, ~ 10 18 seconds. In most cases it has to be even larger than about 
10 26 seconds, in order that the decay doesn't cause fluxes of cosmic positrons, antiprotons and y 
larger than the observed ones. It turns out nevertheless that in most DM models this astonishing 
fact is not explained but rather assumed by hand, typically by assuming a discrete Z2 symmetry. 
That, of course, does not mean at all that these models couldn't be correct, but certainly that their 
completeness is questionable. The origin of this stability is most probably preponderant in deter- 
mining the structure of the DM interactions. Besides providing, for example, an explanation for the 
DM relic density, it would be highly desirable that DM models also provide a fundamental expla- 
nation for this stability. Here, by fundamental, we mean from a dynamical reason, resulting from a 
gauge symmetry and the particle field content of the model. A strong motivation to look for such 
stabilization mechanism is the situation of the Standard Model (SM). In the SM there is a series of 
stable particles whose stability results in all cases from the gauge symmetries and particle content, 
rather than from an ad-hoc symmetry. The photon is stable because it is the massless gauge boson 
of the exact electromagnetic U(\)qed gauge symmetry. The e~ is stable because it is the lightest 
particle charged under this gauge group. The lightest neutrino is stable because of Lorentz invari- 
ance, since it is the lightest fermion. Finally the proton is stable because of conservation of baryon 
number, which is not imposed by hand, but results accidentally from the SM gauge symmetries 
and the gauge charges assigned to the SM particles. In the following we review briefly the various 
simple stabilization mechanisms which can be invoked. We stress that the stabilization mechanism 
has indeed a preponderant influence on the structure of the model and therefore on the associated 
phenomenology. We will distinguish 2 kinds of mechanisms, the low energy ones, which similarly 
to the SM allow to understand the particle stability directly at low energy, and the high energy ones 
which rely on the existence of some explicit UV physics, typically on the existence of a high energy 
gauge group. The latter possibility is less interesting phenomenologically because it relies on an 
assumption difficult to probe experimentally, but is certainly theoretically attractive too, especially 
if it is connected to grand-unification. 

2. The MSSM case 

In the minimal supersymmetric standard model, the lightest supersymmetric particle (LSP), 
in general assumed to be the lightest neutralino, is stable due to the R-parity symmetry. R-parity 
has not been assumed originally to get a stable DM particle but to avoid proton decay. It is an 
attractive feature of the MSSM that DM stability results from a symmetry which has been assumed 
for an other reason. However, yet, R-parity is assumed totally by hand in the MSSM and is not 
explained. Moreover, to avoid proton decay, it would be sufficient to assume the conservation of 
baryon number or lepton number. This would allow terms which do break R-parity and hence 
would destabilize the LSR 

There is nevertheless a way to justify R-parity from a gauge principle. R-parity which is a 
symmetry distinguishing partners and super-partners, is connected to the superfield R-symmetry, 
under which all quark and leptons superfields are odd and the Higgs superfields are even. For any 
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particle of spin s, R p = R s - (— \) 2s . This symmetry prevents any R-parity breaking interactions. 
For any superfield of baryon number 5 and lepton number L, R s turns out to be nothing but R s = 
(— 1) 3 ( B ~ L ). This means that R symmetry is the discrete Z2 subgroup of any gauge U(\)b-l group 
we could have on top of the SM gauge groups. Consequently it is is also a subgroup of the grand- 
unified 50(10) group. Therefore 5-symmetry could be the discrete remnant subgroup of these 
gauge groups. If all scalar fields, whose vacuum expectation values break these groups, are even 
under 5 — L, i.e. have even R s , /^-symmetry will remain automatically unbroken [1, 2], since they 
will leave all interactions even. This is the well-known UV justification for R-symmetry. The 
phenomenology of these models is closely related to the fact that, in order that the mechanism 
works, all scalar representations breaking the various gauge symmetries have to be taken even. 
Natural in these theories are possible new states at the accessible TeV scale, b-x unification, or 
scalar triplet seesaw as dominant origin of neutrino masses [2, 3]. 

3. Non-supersymmetric U(1)b-l justification of DM stability 

During a long time the above GUT justification of DM stability was somehow believed to be 
a privilege of supersymmetric theories. It is interesting to stress that this possible DM stability 
explanation is independent of supersymmetry. The same mechanism can be invoked in the non- 
supersymmetric context. This is a possibility which has been considered recently in Refs. [4, 5, 
9]. In the SM all quarks and fermions have odd 3(5 — L) and in particular in 50(10) an entire 
generation of them can be put in a single 16 v representation which is 3(5 — L) odd. The Higgs 
fields on the other hand are even, and can be put in the 10// representation of 50(10), which 
is even. As a result if all fields breaking U{\)b-l an d 50(10) are 3(5 — L) even (i.e. among 
10//, 45//, 54// , 120//, 126//, 210//, ...), a Z^ 1 symmetry will remain unbroken. This leads to 2 
possible classes of models for DM. Either DM is scalar in case to be stable it has to belong to an 
odd representation of 3(5 — L), because all other scalar particle combinations it could couple to, 
to decay, are even. The 5O(10) odd representations it could belong to are 16//, 144//, 560//, as 
showed in Refs. [4]-[7]. Or DM could be a fermion [9], if it belongs to an even representation, 
among 10^,45^,54^,12011/, 126^,210^, .. .because all other fermion combination it could couple 
to are odd. This leads to two minimal models: either it belongs to a scalar 16// representation, in 
case it is necessarily a combination of the electroweak singlet and the neutral component of the 
electroweak doublet in this representation. These are the only 2 neutral and colorless particles in 
this multiplet, and we know DM has to be neutral and colorless. Or it belongs to the fermion 45^ 
(or very similarly 54^) representation, in case it has to be the neutral component of the fermionic 
electroweak triplet in this representation. 

The phenomenology of the scalar candidate has been considered at length in Refs. [4, 5, 6, 7]. 
The scalar DM particle stands within the range ~ 70 GeV-few TeV, the lower bound coming from 
the LEP constraint on the search for charged Higgs particles. The direct detection cross section rate 
lies typically between 10 cm 2 and 10 cm 2 , therefore in principle detectable in most cases, 
see Fig. 3 of [5]. The next-to-lightest dark scalar as well as the charged scalar DM partner in the 
doublet, are predicted to be long-lived in various regimes, providing experimental signatures in the 
form of displaced vertices with missing (DM) energy. Electroweak symmetry breaking can occur 
radiatively in this setup from the DM-Higgs couplings. 
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The fermion case [9] is in our opinion interesting in many respects. First of all it is quite 
predictive. As the DM is made of the neutral component of the hyperchargeless fermion triplet, 
the only interactions it can have are the gauge interactions (unlike a scalar which can always have 
unknown quartic scalar interactions). It is therefore a pure "WIMP" candidate whose annihilation 
cross section (to gauge bosons or to a fermion pair) is completely fixed by the gauge interactions 
and DM mass. As a result the observed DM relic density 0.091 < Q.DMh 2 < 0.129 [8] can be 
obtained only for a fixed value of the DM mass which turns out to be 2.7 ±0.15 TeV [19, 22, 
20]. Recently, taking into account the effect the DM kinetic decoupling has on the Sommerfeld 
enhancement of the DM annihilation, the right relic density has also been found for nijjM — 4 TeV, 
where a Sommerfeld resonance stands [10]. Secondly since we justify the DM stability by grand- 
unification we must find a way in this model to get unification of the SM gauge couplings at the 
GUT scale, M GUT ~ 10 16 GeV. It is remarkable that the adjunction to the SM of the fermion DM 
triplet pushes the unification scale (of the hypercharge and weak couplings) from ~ 10 12 GeV scale 
to the GUT scale (which is sufficiently high to avoid proton decay problems). In other words if 
DM is a WIMP, which is the most straightforward possibility, we know it has to be around the 
electroweak scale or below. It is therefore fully motivated to believe, on the basis of the DM 
experimental fact, that there should be at least one extra particle on top of the SM ones around the 
electroweak scale. It turns out that, adding just this DM WIMP, one can have unification of the 
gauge couplings. 1 In other words this model succeeds in realizing in a particularly minimal way 
the features split-supersymmetry was mostly made for (i.e. DM and gauge unification), without 
needing to assume any supersymmetry. Of course this model, as split supersymmetry, doesn't 
address the hierarchy problem. However to have an extra fermion around the electroweak scale 
doesn't appear as unnatural in the sense that it is protected by a t/(l) global chiral symmetry in the 
limit of vanishing DM mass. 

The direct detection rate of a 2.75 TeV fermion triplet has been determined in Ref. [19]. 
It is of order 10~ 45 cm 2 , hence in principle reachable. As for indirect detection, a promising 
possibility is to observe, in atmospheric Cherenkov experiment looking at the galactic center [11], 
monochromatic y lines from its annihilation. The corresponding rates for a 4 TeV candidate are 
given in Ref. [10], where the production of positrons has been found to be largely boosted by the 
Sommerfeld effect, as required to reproduce the e + excess observed by the Pamela satellite. 

4. Exact gauge symmetry setup 

In the following we will review 5 setups to justify the DM stability in a way which can be 
understood directly at low energy, i.e. without needing to assume a high energy gauge symmetry. 
Three of them, section 6, rely on the existence of an accidental symmetry. The fourth one, section 

5, relies on the existence of a remnant symmetry of a gauge group broken at low energy. The fifth 
one, in this section, relies on the existence of an exact low energy gauge symmetry. 

In the SM the e~ is stable because it is the lightest charged particle under the unbroken electro- 
magnetic gauge group, U{Y)qed- Similarly if there were, on top of the SM gauge groups, an extra 

'Of course, to achieve unification of the 3 gauge couplings, one also needs as usual to add colorful particles, in 
order to increase the value of the QCD coupling at the GUT scale. This is provided by the setup too, as the 45 (or 54) 
fermion representation contains a color octet which can do the job [9]. 
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Figure 1: Values of a' needed to get the observed relic density, as a function of moM ( m GeV), for m/, = 120 
GeV and £=0.1, 0.5, 1. 

unbroken U (1)qed'> the lightest charged particle under it would be stable. Actually it is impressive 
that the adjunction to the SM of the simplest gauge structure one could think of, that is to say QED 
for a single fermion, 



with = dp — ie'Ap, leads to a viable DM candidate! This QED' model, on top of the SM 
particles, contains an extra massless / and an extra electron, e' (singlet of the SM gauge groups), 
and nothing more! It has been studied in Ref. [12], and in Refs. [13, 14] (in the context of a larger 
framework consisting of a copy of the MSSM in the hidden sector, with the DM in the form of 
a stau). This structure also appears in the context of the mirror models, although in these models 
DM is in general dominated by the mirror baryons rather than by the mirror electron. Unless the 
QED' fine structure constant a' is tiny, DM is in thermal equilibrium with the /. The DM sector, 
however, as considered in Refs. [12, 13, 14], is not in thermal equilibrium at all with the SM 
visible one, since no connector between both sectors are basically considered. The relic density is 
determined by the usual freeze out of its pure hidden sector annihilation cross section e'e' — > //. 
It depends only on 3 parameters: m e <, a! and the temperature of the hidden sector relative to the SM 
one, t, = T' I J '. The S-wave dominated cross section is G ann ih • v ~ net' 2 / (2m| JAf ). Analytically, for 
the freeze-out temperature Tf we obtain: 



(4.1) 



x f = £,-ln[0m%-Z, 5ll O annih vM pl m D M—^] 



-Z - • ln{$ ■ Zrc[0.038 • $ 5/2 o wmih vM pl m DM ^_}} 



(4.2) 



with Xf = mDAt/Tf which gives 
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Here it has been assumed that the present value of ^ 3 is small, so that ^(today);^ 3.91 
(which implicitly means that the Hubble expansion rate H is mostly determined by the visible 
sector particle species). In Fig. 1 we display the values of m e i and a' which lead to a relic den- 
sity within the WMAP range. They can be understood in the following way. Given the non- 
relativistic Boltzmann suppressed behavior of the DM equilibrium number density around freeze- 
out, udm{T') °= (T') 3 / z e~ mDM / T ' , j-jjg va i ue f niDM/T'f = Xf/E, doesn't change much with even 
though the smaller £ , the larger the Hubble constant at T' ~ Tl (when T ann \h ~ H). The dependence 
is only logarithmic. But if one considers a smaller £ the number of DM particles left at freeze- 
out increases as H ~ <^ -2 (since r anni 7,(7^) ~ o ann ihvriDM{Tf))- On the other hand the entropy at 
freeze-out, which is proportional to Tf, increases as <i;~ 3 . Therefore Q-dm = ™dm«om° V / Proday — 
niDM{n DM {T'j-)/s{T'j)).{sToday/PToday) is proportional to E,, see Eq. (4.3). For smaller £ this must 
be compensated by taking a smaller c ann jh, as shown in Fig. 1, so that the freeze-out occurs earlier. 
Note that these results agree with the one of Ref. [13] up to the power 5/2 instead of 3/2 in Eq. (4.2) 
(which is freeze-out prescription dependent and has a moderate numerical effect). 2 

The QED' model is a typical example of model where the mechanism at the origin of DM sta- 
bility implies an all specific phenomenology. The presence of an extra unbroken gauge group 
implies extra radiation in the Universe, i.e. the presence of a long range force which acts on 
the charged DM particles. This has an all series of cosmological and astrophysical implications 
[12, 13, 14], such as a) e'-e' boundstate formation which can change the DM annihilation rate 
at various epochs, b) delaying of the DM kinetic decoupling from yV — > yV Compton scatter- 
ing which can delay the formation of small scale structure, c) DM self interactions induced by 
e'e' — > e'e' Rutherford scattering which may affect the dynamics of the Bullet Cluster as well as 
the morphology of DM galactic haloes. The main constraint turns out to be the latter: more colli- 
sions affects the ellipticity of the DM halo, tending to the formation of a more isothermal profile, 
that is to say to the formation of a core, a region of constant DM density around the galactic 
center. It gives an upper bound on a' which decreases if mdm decreases. The upper bound is 
a' < 1(T 8 , lfr 5 , 1(T 2 for m DM = 100 MeV, 10 GeV, 1 TeV respectively. Combined with the relic 
density constraint above it gives the upper bound £ < 10 5 , 10~ 3 , 10 _1 respectively. In compari- 
son, the BBN constraint on extra relativistic degrees of freedom at Tbbn ~ 1 MeV is very mild. If in 
the hidden sector all particles have a mass exceeding this temperature, except the /, the constraint 
is 0.9. 

Given these small values of E, , the connectors between this QED' structure and the SM visible 
sectors must be small, otherwise both sectors would thermalize, giving E, ~ 1 . For the fermionic 
case of Eq. (4.1) the only possible (i.e. gauge invariant) term is a kinetic mixing term J2? B 
kF^ v F^ v , with F£ v the hypercharge gauge field. One can determine that K has to be smaller than 
10 9 — 10 1 1 to avoid a too large transfer of energy from the visible to the hidden sector (compared 
to the ellipticity bound on above). This unfortunately prevents us to have access to the hidden 
sector and DM at colliders. 

In the case of a scalar DM candidate, i.e. scalar electrodynamics for a single hidden sector <p' 
field, the phenomenology is very similar, apart for one important difference: this DM scalar can 
directly interact with the SM through the Higgs portal Xmty'ty'^WH interaction. This interaction 

2 These results have been obtained in collaboration with Xiaoyong Chu, to appear in a longer publication. 
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has nevertheless to be small for the same reason as for K. Unfortunately, simplicity doesn't always 
mean testability at colliders. 

5. Remnant symmetry of a gauge group broken at low energy 

In the above one has seen that DM stability could result from the existence of a remnant global 
symmetry subgroup of a high energy gauge group. This requires that the gauge group is broken by 
scalar fields which are even under the remnant group. This, actually, could hold at low energy too, 
either by considering that U (1 ) b-l is broken at low energy, or independently of B — L. The simplest 
possibility in the latter case is to consider an extra U(l)' gauge group with 2 particles charged under 
it, a scalar (j)', which breaks it spontaneously, and a fermion y'. The relevant lagrangian is 

JSf = J^m + ^ V v + Yiip - m ¥/ )xi/ + D^'D^ n 

-1%$' <j>' t -V(ff t ) 2 -^>' t ^ t +KF^ V , (5.1) 

with H the SM Higgs doublet and with fli, < 0, so that (j)' acquires a vev, (0') = v^. After symmetry 
breaking the model contains the DM fermion, a real scalar, v\ , and a massive Z' gauge boson. Since 
0' has a vev, v\ is obviously not stable but y' remains stable even though £7(1)' is broken, simply 
because all interactions involve it in pairs. Such a model has been considered in Ref. [15], in the 
so-called "secluded" DM framework. Since the gauge group is broken, the U(l)' gauge boson 
is massive. Therefore all the constraints due to extra radiation in the QED' model above are not 
relevant anymore. The interactions connecting the hidden and visible sectors, kinetic mixing and 
Higgs portal, can be much larger than for the unbroken QED' model above. For instance K can be 
as large as few 10~ 3 . The relic density can proceed from the purely hidden sector process involving 
the Z', in case it is "secluded" or, if the connectors are large, can result from the annihilation induced 
by these connectors. Interestingly, through the exchange of a Z', the kinetic mixing interaction fc 
induces an interaction between a DM pair and a SM fermion pair, y'y— > Z'— > y,Z — > /sm/sm- 
If the DM mass is of order 10 GeV, one can get [16] from this single diagram both the observed relic 
density and a direct detection rate of order the ones indicated by the Dama and Cogent experiment 
[17, 18]. This requires nevertheless a rather precise correlation between mdm and m z /: the mass of 
the Z' must be within a range close to twice the DM mass value, in order that the DM annihilation 
is largely enhanced by a Z' resonance. 

6. Accidental stability setups 

As said above there exists already in the SM an example of particle which is accidentally 
stable, the proton. There is no reason why DM could not be stable in a similar way. Actually one 
would expect in this case an interesting property to hold: decay induced fluxes of cosmic rays of 
order the ones observed. If the stability is accidental, there is in particular no reason why DM 
should be absolutely stable, unlike in the DM models considered above. Even if the pure low 
energy renormalizable interactions with the SM particles do not allow DM to decay, in this case 
there is no symmetry protecting DM from being destabilized by an interaction involving a higher 
energy particle beyond the SM. On the contrary, if from the exchange of such a heavy particle, a 



7 



On the stability of particle dark matter 



Thomas Hambye 



dimension 5 interaction (suppressed by one power of the heavy particle scale A), is induced, it will 
induce a far too short lifetime. The expected DM lifetime in this case is Xdm ~ 87rA 2 /ra 3 where 
m can only be a combination of the DM mass and SM particle masses. Considering a DM mass 
below few TeV leads to a lifetime shorter than the age of the Universe, %u, unless A > 10 25 GeV, 
i.e. far above the Planck scale. An even larger scale ~ 10 29 GeV is needed, if one demands that the 
lifetime is larger than ~ 10 sec, what is mandatory for most models in order that the DM decay 
does not induce a flux of cosmic rays larger than the ones observed. One therefore must make sure 
that no dimension five operator between the DM candidate field and SM fields is allowed by the SM 
(and hidden sector) gauge symmetries. But if, on the contrary, they allow dimension 6 operators it 
turns out that the lifetime, Zdm ~ 87rA 4 /m 5 , is of order 10 26 sec if A ~ 10 14 GeV, which is close to 
the GUT scale! In other words an accidentally stable DM candidate which can be destabilized by a 
dimension 6 GUT scale induced interaction leads to a flux of cosmic rays of the order of the ones 
observed, and therefore potentially to a rich phenomenology. 

6.1 Minimal Dark Matter 

Probably the simplest setup which leads to an accidentally stable DM candidate is the Minimal 
DM one [19, 20, 21]. This framework doesn't require any extra gauge groups on top of the SM ones 
or the assumption of a global group. It requires, instead, high SU (2)i multiplets. It stems from 
the fact that if one considers a large enough fermion or scalar SU (2)l multiplet, SU {2)i invariance 
prevents it to have any "fast" renormalizable or dimension 5 interactions. For instance, to have 
a renormalizable interaction involving a fermion SU {2)i quadruplet, one needs it to couple to a 
quadruplet of dimension 5/2 made of SM fields. This can only involve a fermion, which in the SM 
belongs at most to a doublet, and the SM Higgs doublet scalar. From 2 doublets one gets at most a 
triplet. Therefore a fermion quadruplet or higher has no renormalizable interactions with SM fields, 
but still can have a dimension 5 interaction (with one lepton doublet and 2 Higgs doublets). But 
a fermion quintuplet or more cannot have a dimension 5 interaction either. Similarly, for a scalar 
candidate, a sextuplet or higher cannot have any of these interactions. 

To this constraint come 2 other ones. First DM must be neutral, hence the multiplet must have a 
neutral component, i.e. an integer hypercharge. Second this neutral component cannot have any tree 
level interactions to the Z boson, otherwise the direct detection rate driven by Z exchange would 
be much larger than allowed experimentally. In other words the electric charge Q must vanishes as 
well as the hypercharge Y and the weak isospin Tt„ which is possible only for an odd multiplet. The 
good candidates are therefore a fermion quintuplet, septuplet,... or a scalar septuplet, nonuplet,... 
Why there would be beyond the SM so large multiplets, with no other TeV scale smaller multiplets 
(which would allow to write again fast destabilizing interactions), is not clear but this is certainly a 
possibility. 

We will not review here the phenomenology of such a candidate. It has been summarized in 
Ref. [21]. Some of its main features are that a) the fermion quintuplet must have a 9.6 TeV mass 
in order to satisfy the relic density constraint, b) the direct detection rate is of order 10 cm 2 , c) 
the indirect detection rate is boosted by two orders of magnitude by the Sommerfeld effect, so that 
the Pamela e + excess can be reproduced if one assumes an extra astrophysical boost of order 50 
(but giving in this case a flux of e + + e~ exceeding the Hess experiment value at energies above 
~ 2 TeV), d) given the high DM mass this doesn't come with a too large flux of p (because the 
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excess of p it gives is above 100 GeV), and e) the expected flux of /produced in the galactic center 
tends to exceed the experimental limits in the case of a NFW halo profile, but a slightly more 
isothermal profile is fine. 

6.2 Hidden vector DM 

Another setup based on an accidental stabilization is the hidden vector DM model [23, 24, 25]. 
This model in many respects is different from the other models above. First, this model shows that 
a non-abelian gauge boson could be a DM particle candidate. Second, its stability results from the 
automatic presence of an accidental symmetry which is non-abelian, unlike all models above which 
are based on a discrete or global U(l) symmetry. Third, DM can be destabilized only by dimension 
6 or higher operators. No need to impose that there are no dimension 5 operators. Fourth, related to 
the fact that DM has spin-1, the dimension 6 operators generically lead to intense monochromatic 
/ray lines directly at tree level, rather than at the usual suppressed one-loop level. As well known 
the observation of such lines would be a real smoking gun for DM, because there is no expected 
relevant astrophysical background for such a signal. 

This model works as follows. Consider a non abelian gauge structure on top of the SM ones. 
The simplest example is to consider a SU (2) one. Consider also a scalar multiplet of it, which 
breaks SU (2) completely, in such a way that all 3 gauge bosons becomes massive (no extra radiation 
in the Universe). The most general Lagrangian for such a structure is 

Sf = J? SM - ±F ,tiV • f; v + (D M 0) t (D» - ;U W# - A^ + - W , (6.1) 

with D^(j) = d^tj) — i^-x • A /fX . It contains a pure gauge term with gauge coupling g^, a kinetic term 
and potential for the scalar and a Higgs portal interaction. The latter is not forbidden by any of 
the visible or hidden sector gauge symmetry. After <p gets a vev, , the 3 massive gauge bosons 
have a mass my = g§v<$, /2 and a single real scalar Higgs boson, r\, remains. The all model can 
be parametrized in terms of 4 parameters: g§, my, m-q and X m . Looking closer at this structure, 
which is similar to the scalar-gauge SM one (but without extra £7(1), i.e. no Weinberg angle), what 
one remarks is that, even if the gauge group has been broken completely, a non-abelian global 
symmetry remains accidentally: a custodial 50(3) group under which the 3 massive gauge bosons 
form a triplet vector and the tj a singlet. Since this symmetry is exact, any non singlet particle 
cannot decay into singlets, hence the massive gauge bosons are stable. Related to the custodial 
symmetry, the 3 gauge bosons have the same mass. The custodial symmetry of the pure gauge- 
scalar hidden sector is not broken either by the Higgs portal interaction, i.e. it can couple to the SM 
without being destabilized. 

The massive vectors in this model can fulfill the various DM constraints. Their relic density 
results from the thermal freeze-out of their annihilation. We can distinguish 2 regimes depending 
on the size of the Higgs portal interaction. If X m is small the relic density results from the pure 
hidden sector annihilation, V/V; — > T]7] and V/V) — > ViJ]. Note that the last annihilation process is a 
2 to 1 DM process, impossible in model based on a usual Z2 symmetry. It is induced by the non- 
abelian gauge trilinear term. If A,„ is instead large, typically larger than 10~ 3 , extra annihilation to 
SM particles can be important or even dominates. Part of these channels comes from the fact that 
X m induces a mixing of the 77 boson with the SM Higgs boson h. Values of m r] and my leading to 
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the right relic density in both regimes cover wide ranges of values. Direct detection rates saturating 
the present experimental bounds as well as signatures at LHC can be easily obtained in the large 
Higgs portal regime. For more details see Refs. [23, 25]. Here instead we would like to discuss 
more the monochromatic y flux it is expected to lead to. 

Since the custodial symmetry is accidental it could be broken, as explained above, by interac- 
tions induced from the exchange of higher energy fields. These interactions, which involve only 
low energy particles as external particles, have obviously to respect the gauge symmetries of both 
sectors. Nicely there exist no dimension 5 operators of this kind, but 4 dimension 6 ones which 
break the custodial symmetry: 

(A) ^D^<j>D^H, (6.2) 

(B) ^D^^H^D^H, (6.3) 

(C) ^D^D^F^ 7 , (6.4) 

(D) ^fF^F^ 7 . (6.5) 

Associated to the spin 1 character of DM, all of them turn out to generically induce, with large 
branching ratios, 2 body decays involving a photon in the final state (from F^ vY or from the covari- 
ant derivatives of the SM Higgs doublet). As explained above, if the underlying high physics scale 
is of order the GUT scale one expects fluxes of cosmic rays of the order of the present sensitivities. 
This is shown in Fig. 2 which, together with the flux of other cosmic rays produced, gives an ex- 
ample of monochromatic lines one can obtain assuming such a scale. Note that, since these lines 
proceed through a decay, their magnitude is proportional to the number of DM particles «dm> rather 
than to rv\ )M as for an annihilation. As a result one doesn't need to look at the "involved" galactic 
center to expect an observable monochromatic line. The flux given in Fig. 2 is the extragalactic 
flux. Such monochromatic lines could be observed by the Fermi satellite if it increases its statistics, 
or could extend its analysis towards higher energies. 

Finally, about this model, note that if, instead of being broken by the vacuum expectation 
value of the scalar field 0', the SU (2)hs gauge structure of Eq. (6.1) confines at low energy (i.e. if 
there is no (j)' vev or if its value would be smaller than the confinement scale) the model leads to 
accidentally stable spin 1 states which are good DM candidates too [24]. In this case the stable 
vectors turn out to be boundstates of the scalar doublet and gauge bosons, of the type 0^D^0. They 
are stable because they also form a triplet of the custodial symmetry, and custodial symmetry is not 
broken by the confinement. Given the strongly interacting interactions it involves, DM particle is 
expected to stand in the multi TeV range, to satisfy the relic density constraint. 

6.3 Weakly interacting stable pions 

Another possibility of accidentally stable DM candidate arises in the context of strongly inter- 
acting theories, as the last case discussed above, but in a different way. DM in this case are scalar 
boundstates made of fermions, similar to pions in the SM [26]. Such a structure can arise if on top 
of the SM there is a new confining structure, a "QCD"' sector, under which there are new fermions 
which are vectorlike, i.e. whose chiral components couple in the same way under the new strong 
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Figure 2: Predictions [25] for m v = 1550 GeV, g$ = 2.1, = 1457 GeV, m n = 1243 GeV, m h = 153 GeV, 
and for the rj — h mixing, sin j3 = 0.25, considering as an example only operator C, with z^m = 1 -6 x 10 27 s 
(A= 1.2 X 10 16 GeV). The upper panels show the positron fraction (left) and the total electron + positron 
flux (right) compared with experimental data. Dashed lines show the adopted astrophysical background, 
solid lines are background + dark matter signal (which overlap the background in this plot). The lower 
left panel shows the gamma-ray signal from dark matter decay, whereas the lower right panel shows the 
pi p-ratio: background (dashed line) and overall flux (solid lines, again identical with background). 



interactions (as quarks for QCD) as well as under weak interactions (unlike quarks in the SM) 

& = &sm ~ \F%F» va + fyW +g 4h b + g w f a J a )Wi (6-6) 

where = l,...,n are the new fermions, t" are the generators of the new strong group and J" 
are the SM generators of SU (2)l- As in QCD, after the associated chiral symmetry gets dynami- 
cally broken, one obtains a hadronic mass spectrum whose lightest states are the chiral symmetry 
Goldstone bosons, the pions, "uf". As QCD too, the QCD' interactions respect isospin and charge 
conjugation, therefore respect the G-parity. G-parity conservation in QCD explains why a p meson 
which has an even G-parity can decay strongly to 2 %, but not to 3 pions, since a pion has odd 
G-parity (and the contrary for the ft) meson which has odd G-parity). In the SM G-parity is nev- 
ertheless broken by the electroweak interactions due to the chiral structure of quarks under them, 
so that the lightest G-odd states, i.e. the pions, are not stable. Here on the contrary since these 
new strongly interacting fermions are assumed to have the same electroweak interactions for left 
and right-handed components, G-parity is not broken and the %' are stable. The lightest %' cannot 
annihilate strongly, since it is the lightest hadron, but do annihilate weakly. Therefore it is a real 
WIMP and has the right relic density if its mass is of order TeV. As a result such a model can be 
probed at the LHC. For the direct detection the cross section is of order 10~ 45 cm 2 , as for typical 
SU(2)l multiplets, see section 6.1. 
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Since G-parity is a purely accidental symmetry these pions could in principle decay from the 
exchange of heavy states. A problem one encounters in this framework is that, having a vectorial 
structure to guaranty the conservation of G-parity, one can write down G-parity violating dimension 
5 operators destabilizing the pions, of the form F^ v \j/a^ v y/A and H^J a H\j/J a \j//A. The only way 
to avoid them appears to assume an extra global symmetry, such as a Peccei-Quinn symmetry. This 
weakens the naturalness of the stability of such candidates, but is certainly a possibility too. 

Before concluding, note that another accidentally stable DM candidate arises in the context 
of the mirror models [27], in the form of mirror nuclei, protected by the conservation of mirror 
baryon number. Note also that in the above we focus our attention on dynamical ways to justify 
the DM stability, because the gauge symmetry principle constitutes the most established and robust 
SM pillar. No global symmetry have still been firmly established. However, it is highly probable 
that behind the quark and lepton flavor puzzle there exist underlying flavor symmetries. Another 
possibility is therefore that the DM stability is closely related to the flavor symmetry structures, for 
instance that the discrete symmetry at the origin of DM stability would be a remnant symmetry of a 
broken flavor group. This possibility has been recently considered [28] on the basis of an A4 global 
flavor symmetry which is dynamically broken to a the Z2 discrete group. 

Finally we limited ourselves to the candidates whose relic density is determined by the thermal 
freeze-out of their annihilation. Of course there also exist DM candidates whose stability results 
from the fact that the allowed renormalizable couplings, which could induce a decay, would be very 
tiny. In particular if the only interaction causing DM decay is the gravitational coupling, one can 
have a stable enough DM candidate. Similarly if one assumes very light right-handed neutrinos, 
of order KeV, and therefore tiny Yukawa couplings (not to induce too large neutrino masses) the 
lightest right-handed neutrino could be stable enough [29]. Another DM candidate is of course 
the axion which if it has a very tiny mass, ^ ( 10 3 eV), can decay only to photons with naturally 
suppressed rates [31]. A KeV majoron is also an option [32]. 



7. Conclusion 

The fact that the DM particle is stable on cosmological time scales is a peculiar property. It 
definitely needs an explanation. On the basis of a gauge principle there are quite a few possible 
origins for this stability. Along these mechanisms DM can be a scalar, a fermion, or even a gauge 
boson (or higher spin object such as the gravitino). The gauge symmetry invoked can be abelian or 
non abelian (both in their confined or unconfined phase). Each mechanism leads to a characteristic 
phenomenology. In particular, 1) the existence of an exact gauge group results in extra radiation 
which might have many astrophysical effects, 2) accidental symmetry mechanisms are expected to 
lead to a rich indirect detection phenomenology from DM decay, such as intense 7 ray lines in the 
hidden vector DM setup, 3) high energy stabilization mechanism, based in particular on 50(10), 
can have DM as the missing piece for unification of gauge couplings, 4) if no extra gauge group 
are assumed the DM must be part of a high SU (2)l multiplet with definite properties. Apart for 
the models where DM lies in the multi-TeV range, and for the models based on extra-radiation, all 
models can lead to specific signatures at the LHC. 
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